Pb adatom activated dissociation of O2 and oxidation of Pb 

surfaces 



Yu Yang/'^ Jia Li,^ Zhirong Liu,^ Gang Zhou,^ Jian Wu,^ Wenhui Duan,^* 
Peng Jiang,^ Jin-Feng Jia,^ Qi-Kun Xue,^'*^ Bing-Lin Gu,^ and S. B. Zhang^ 
^LCP, Institute of Applied Physics and Computational Mathematics, 
P.O. Box 8009, Beijing 100088, People's Republic of China 
'^Department of Physics, Tsinghua University, 
Beijing 100084, People's Republic of China 
^College of Chemistry and Molecular Engineering, 
Peking University, Beijing 100871, People's Republic of China 
"^Institute of Physics, Chinese Academy of Sciences, 
Beijing 100080, People's Republic of China and 
^Department of Physics, Applied Physics, and Astronomy, 
Rensselaer Polytechnic Institute, Troy, NY 12180, USA 
(Dated: August 24, 2009) 

Abstract 

We investigate the dissociation of O2 on Pb(lll) surface using first-principles calculations. It is 
found that in a practical high vacuum environment, the adsorption of molecular O2 takes place on 
clean Pb surfaces only at low temperatures such as 100 K, but the O2 easily desorbs at (elevated) 
room temperatures. It is further found that the Pb adatoms enhance the molecular adsorption 
and activate the adsorbed O2 to dissociate during subsequent room temperature annealing. Our 
theory explains the observation of a two-step oxidation process on the Pb surfaces by the unique 
role of Pb adatoms. 
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It is of fundamental importance to understand the adsorption, dissociation, and the corre- 
sponding kinetic processes of diatomic molecules on materials surfaces As a prototypical 
example, the adsorption and dissociation of O2 gas on metal surfaces have attracted con- 
siderable attention because the resulting metal oxides have been widely used as catalysts, 
sensors, dielectrics, and corrosion inhibitors In this regard, lead (Pb) (111) surface is a 
good model system, because atomically-fiat Pb(lll) terraces have been fabricated on both 
metal and semiconductor substrates such as Ru (0001) [3] and Si(lll) Rich experimen- 
tal and theoretical literatures are now available for such surfaces sl- It was found that O2 
molecules adsorb on Pb(lll) at temperatures as low as 100 K. Subsequent room temperature 
annealing will not result in the desorption of O2 but the oxidization of the Pb surface Gj. In 
contrast, the Pb(lll) surface is remarkably resistant to oxidation at room temperature in 



the experimental environment 



,|7|]. From simple energetic argument, it is very difficult to 
reconcile these experimental observations, namely, how can the molecularly adsorbed O2 at 
low temperature does not desorb at (elevated) room temperature annealing, while the same 
O2 molecules would hardly adsorb on to the same surface at the same temperature. 

In this paper, we explore the mechanisms for the puzzling adsorption and dissociation of 
O2 on Pb(lll) surfaces, based on first-principles total-energy calculations and thermody- 
namic analysis. We show that one cannot determine the state of O2 at different temperatures 
solely from the calculated adsorption energy but rather one has to rely on the Gibbs free en- 
ergy that explicitly includes the contributions from gas-phase oxygen in terms of its atomic 
chemical potential. Furthermore, we show that Pb adatoms, a common surface defect on the 
Pb terraces, play a key role in retaining the O2 adsorbed at low temperatures and assisting 
its dissociation once the temperature is increased. Our results are in good agreement with 
experiments. 

Our calculations are based on the density functional theory as implemented in the Vi- 
enna Ab-initio Simulation Package (VASP) a]. The Perdew-Wang-91 (PW91}j9| generalized 
gradient approximation and the projector- augmented wave (PAW) potential[10| are used to 
describe the exchange-correlation energy and the electron-ion interaction, respectively. The 
nudged elastic band (NEB) method ll|] is used to find the minimum energy path and the 
transition state for O2 dissociation. We confirm that the transition states have only one 
imaginary frequency. The Pb(lll) surface is modeled by a 2 x 2 surface slab with 4 atomic 
layers and a vacuum layer of 15 A thickness. The integration over the Brillouin zone is 
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carried out by using the Monkhorst-Pack scheme [12] with a 7 x 7 x 1 grid. The cutoff en- 
ergy for the plane wave expansion is 400 eV. All the atoms, except for those at the bottom 
layer, are relaxed until the forces are less than 0.01 eV/A. The calculated lattice constant 
for bulk Pb and the bond length for isolated O2 are 5.03 A and 1.24 A, respectively, in 



good agreement with the experimental values of 4.95 A 13] and 1.21 A[lJ]. To calculate the 
adsorption energy, we use the revised Perdew-Burke-Ernzerhof (RPBE) functional for the 
exchange- correlation energy at the already optimized geometries and the self-consistent elec- 
tron densities given by PW91. This scheme is known to yield better results when compared 



with experiments 151 ]. 



We have previously studied [16(] in-depth the molecular O2 adsorption on Pb(lll), which 
establishes the necessary initial state for the O2 dissociation. Figure 1 shows the calculated 
dissociation path where the initial state is the molecular adsorption (MA) state and the 
final state is the atomic adsorption (AA) state. The MA state has no energy barrier for 
adsorption [l^, whereas the AA state is 2.29 eV lower in energy than the MA state. In the 
AA state, the two oxygen atoms reside at two adjacent surface hollow sites. The calculated 
energy barrier going from MA to AA is 0.32 eV. At the transition state (TSl), the two 
oxygen atoms reside on two adjacent surface bridge sites separated by 1.80 A. 



Based on the van't Hoff-Arrhenius law in the harmonic approximation 17|], the rate 
coefficient k for the dissociation can be expressed as 



T-r3Af ,,MA 

^ = nfe^ X = - X (1) 

ffi=i 

where uf^^ and uj^^ are the frequencies of the eigenmodes at molecular adsorption and tran- 
sition states, respectively, z/ is the attempt frequency (for O2, the calculated u = 1.33x10^^ 
s~^), AE is the energy difference between the initial and transition states, and ks, T, and 
are the Boltzmann constant, temperature, and number of atoms, respectively. Using the 
barrier in Fig. 1 {AE = 0.32 eV), we obtain the rate coefficients for O2 dissociation: namely, 
~10'' s^^ at 300 K and ~10~^ s~^ at 100 K. Hence, molecular adsorption is stable at 100 K 
but the adsorbed O2 easily dissociates at 300 K. This conclusion, however, contradicts with 
experiments showing that Pb(lll) surfaces are resistant to oxidation at room temperature 

m. 

The origin of this discrepancy is the crucial role of thermodynamic effect in real experi- 
ments at high vacuum and variable temperature: the Gibbs free energy G{T,p), instead of 
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FIG. 1: (Color online). O2 energies on clean Pb (111) surface. Here, MA and AA stand for 
molecular and atomic adsorptions, respectively, and TSl stands for transition state between the 
two. The corresponding atomic structures are given in the top panel of the figure in which grey 
balls are Pb and red balls are O. The lower inset shows the calculated Gibbs free energies for a free 
O2 at T = (red dotted line), 100 (blue dotted line), and 300 K (green dotted line), respectively. 
The oxygen partial pressure in the calculation is 2 x 10^'' Torr. 



the total energy, should be used to determine the stability of the adsorption structures. In 
other words, we need to include the atomic chemical potential of oxygen, fio- Within the 
ideal-gas law, /io can be calculated by lo], Q 



l^o{T,p) = iJo{T,po) + -kBT\n{p/po) 



(2) 



where p is the pressure and Pq = 1 atm. One can find ^o{T,pq) as a function of temperature 



in Ref. 



19|: at 100 K, fio{T,Po) = -0.08 eV and at 300 K, yUo(T,po) = -0.27 eV. For the 



experimental pressure of p ~ 10^^ Torr 6|, [Zl], the calculated fio at 300 K is approximately 
—0.55 eV. When taking into account the effect of the /iq, the energy of the gas-phase O2 is 
1.11 — 0.38 = 0.73 eV lower than that of MA, according to Fig. 1 and its inset. In other 
words, molecular adsorption will not take place at room temperature. The rate coefficient 
at 300 K for a gas-phase O2 to dissociate with an estimated energy barrier of 1.05 eV is 
quite small (about 10^^ s^^). So neither molecular nor atomic adsorption of oxygen is likely 
to occur at this temperature. These results explain the observations that Pb(lll) surfaces 

riQ 

resist oxidation at room temperature p, 17| ■ 
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From the thermodynamic viewpoint the discrepancy at room temperature is resolved, but 
a new problem arises at low temperatures. According to the free energy calculation, at T = 
100 K and p ~ 10^'' Torr, the molecular adsorption state of —0.38 eV is only 0.02 eV lower 
than that of a free O2 (as shown in the inset of Fig. 1). Therefore, at this temperature the 
adsorbed O2 will desorb easily. With AE = 0.32 eV in Fig. 1, the rate coefficient at 100 K of 
~10~''s~^ is very small. Thus the adsorbed O2 has practically no chance to dissociate to stay 
on the surface before it desorbs. This result, however, contradicts with the observation that 
O2 adsorbs on the Pb surfaces at around 100 K and oxidizes the surface upon subsequent 
room temperature annealing Gj. 

To explain the experiment, we note that surface impurities can act as nucleation centers 
for PbO grains from which the grains grow rap idly and auto-catalytically Pb adatoms 



are commonly seen on epitaxial Pb surfaces 2l| and readily diffuse at temperatures as low as 
100 K 2^. Therefore, Pb adatoms could be the missing link between theory and experiment. 



For the adatoms to affect the oxidation, two scenarios should be considered: i) a diffusing 
Pb adatom gets to the vicinity of an adsorbed O2 and the two then interact or ii) an O2 
directly adsorbs near a Pb adatom. Because the Pb adatoms are fast diffusers that very 
effectively transport Pb atoms from step edges to other parts of the terraces, we feel that 
the first scenario could be more relevant to the experiment. 

To study the Pb-02 interaction, we place a Pb adatom at a distance reasonably away 
from the O2 yet close enough for the two to attract each other, for example, at the second 
nearest neighbor bridge sites around the 02- Our previous study showed that an O2 on 
Pb(lll) has two nearly degenerate states IGj, and correspondingly, six different local Pb 
configurations. As a result of the atomic relaxation, they all converge to one of the two final 
states with shorter Pb-02 separations. Figure 2(a) shows the lower-energy configuration, 
denoted as the non- dissociative adsorption (NDA) state. In the NDA state, the 0-0 bond 
length increases from 1.24 A for isolated O2 to 1.50 A. This suggests that the 0-0 bond, 
although not broken, is significantly weakened with a 21 % elongation. The calculated Pb- 
02 binding energy is 0.70 eV, indicating that the adatoms can substantially increase the O2 
adsorption. 

To determine the final state of O2 dissociation in the presence of Pb adatom, we started 
from the NDA state and increased the distance between the two O atoms within a plane 
parallel to the Pb surface plane until there is no restoring force between them, and then 
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FIG. 2: (Color online), (a) O2 energies on Pb (111) surface in the presence of a Pb adatom. 
MA stands for the molecular O2 adsorption far away from any Pb adatom, NDA and DA stand 
for nondissociative and dissociative adsorptions, respectively, around the Pb adatom, and TS2 is 
the transition state between the two. Inset has the same meaning as the inset in Fig. 1, but for 
configurations shown in Fig. 2. (b) and (c) Contour plots of the electron localization function (see 
text) for NDA and DA, respectively, in a plane containing the Pb adatom and two oxygen atoms. 
Grey, blue, and red balls stand for surface Pb, Pb adatom, and O atom, respectively. 

performed structural optimization. Figure 2(a) (far right) shows the lowest-energy final 
state, denoted as the dissociative adsorption (DA) state. In the DA state, the distance 
between two O atoms is 2.94 A, which is 137% longer than the normal 0-0 bond length 
of 1.24 A. Electron localization function (ELF) [23! calculation confirms the dissociation. 
Figures 2(b) and 2(c) show the ELFs for the NDA and DA states. One can see clearly 
weak) covalent bonding in Fig. 2(b), but not at all in Fig. 2(c). Bader topological analysis 



2J] further reveals that the two O atoms in the DA state gain 1.99 and 1.80 electrons. 



respectively, from the neighboring Pb atoms. The interactions between the O's and Pb in 
the DA state are thus ionic. A slight asymmetry in the charge transfer is because the two 
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oxygen atoms are on inequivalent surface sites. 

Figure 2(a) shows the dissociation path for O2 in the presence of Pb adatom. It can 
be seen that in this case, the dissociation barrier of 0.28 eV is much lower than the energy 
difference of 0.70 eV between the NDA and MA states. Note that at 300 K, the MA state will 
evolve to the desorption state (discussed earlier). Hence, an adsorbed O2 is more likely to 
dissociate than to desorb from the surface in this case. The calculated attempt frequency for 
Fig. 2 is u = 2.97x 10^^ s^^ and the corresponding rate coefficient at 300 K from NDA to the 
transition state TS2 is ~10'^ s~^). Therefore, activated by the Pb adatoms, the adsorbed O2 
can easily dissociate during a room temperature annealing. These results bring our theory 
into complete qualitative agreement with experiments [g], 

In summary, we have developed a comprehensive first-principles theory for the dissocia- 
tion of O2 on Pb(lll) surface. It is found that one needs to include oxygen atomic chemical 
potential given by the experimental conditions (i.e., partial pressure and ambient tempera- 
ture) to qualitatively account for the observation that O2 can adsorb on clean Pb surfaces 
at 100 K but cannot at room temperature. In addition, Pb adatoms, commonly seen on 
Pb surfaces, can result in a substantial enhancement of the O2 adsorption. As a result, O2 
adsorbed at low temperatures such as at 100 K is more likely to dissociate than desorbing 
upon subsequent room temperature annealing. These results explain the recent puzzling 
experimental observations as well as demonstrating the important role of surface defects in 
activating the oxidation of surfaces. 

This work was supported by the Ministry of Science and Technology of China (Grant 
Nos. 2006CB605105 and 2006CB0L0601), and the National Natural Science Foundation of 
China. 



[1] G. R. Darling and S. Holloway, Rep. Prog. Phys. 58, 1595 (1995). 

[2] V. E. Henrich and P. A. Cox, The Surface Science of Metal Oxides (Cambridge University 

Press, Cambridge, 1994). 
[3] K. Thiirmer, J. E. Reutt-Robey, E. D. Williams, M. Uwaha, A. Emundts, and H. P. Bonzel, 

Phys. Rev. Lett. 87, 186102 (2001). 
[4] M. H. Upton, C. M. Wei, M. Y. Chou, T. Miller, and T. C. Chiang, Phys. Rev. Lett. 93, 



7 



026802 (2004); H. Y. Lin, Y. P. Chiu, L. W. Huang, Y. W. Chen, T. Y. Fu, C. S. Chang, and 
T. T. Tsong, Phys. Rev. Lett. 94, 136101 (2005). 
[5] A. Menzel, M. Kammler, E. H. Conrad, V. Yeh, M. Hupalo and M. C. Tringidcs, Phys. Rev. 
B 67, 165314 (2003); P. J. Feibelman, Phys. Rev. B 62, 17020 (2000); C. M. Wei and M. Y. 
Chou, Phys. Rev. B 66, 233408 (2002); M. Milun, P. Pervan, and D. P. Woodruff, Rep. Prog. 
Phys. 65, 99 (2002). 

[6] X. C. Ma, P. Jiang, Y. Qi, J. F. Jia, Y. Yang, W. H. Duan, W. X. Li, X. H. Bao, S. B. Zhang, 
and Q. K. Xue, Proc. Natl. Acad. Sci. USA 104, 9204 (2007). 

[7] K. Thiirmer, E. Wihiams, and J. Reutt-Robey, Science 297, 2033 (2002). 

[8] G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 (1996) and references therein. 

[9] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
[10] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
[11] G. Mills, H. Jonsson and G. K. Schenter, Surf. Sci. 324, 305 (1995). 
[12] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 
[13] R. W. G. Wyckoff, Crystal Structures (Wiley-Interscience, New York, 1965). 
[14] K. P. Huber and G. Herzberg, Constants of Diatomic Molecules (Van Nostrand, New York, 
1979). 

[15] B. Hammer, L. B. Hansen, and J. K. N0rskov, Phys. Rev. B 59, 7413 (1999). The use of 
the non-self-consistent electron densities introduces a negligible error due to the variational 
principle. 

[16] Y. Yang, G. Zhou, J. Wu, W. H. Duan, Q. K. Xue, B. L. Gu, P. Jiang, X. C. Ma, and S. B. 
Zhang, J. Chem. Phys. 128, 164705 (2008). 

[17] G. H. Vineyard, J. Phys. Chem. Solids 3, 121 (1957). 

[18] T. Vegge, Phys. Rev. B 70, 035412 (2004); H. Lee, J. Li, G. Zhou, W. H. Duan, G. Kim, and 

J. Ihm, Phys. Rev. B 77, 235101 (2008). 
[19] K. Renter and M. Scheffler, Phys. Rev. B 65, 035406 (2001). 

[20] Y. F. Zhao, Y. H. Kim, A. C. Dillon, M. J. Heben, and S. B. Zhang, Phys. Rev. Lett. 94, 

155504 (2005). 

[21] S. H. Chang, W. B. Su, W. B. Jian, C. S. Chang, L. J. Chen, and T. T. Tsong, Phys. Rev. B 
65, 245401 (2002); C. S. Jiang, S. C. Li, H. B. Yu, D. Eom, X. D. Wang, P. Ebert, J. F. Jia, 
Q. K. Xue, and C. K. Shih, Phys. Rev. Lett. 92, 106104 (2004). 

8 



[22] T.-L. Chan, C. Z. Wang, M. Hupalo, M. C. Tringides, and K. M. Ho, Phys. Rev. Lett. 96, 
226102 (2006). 

[23] Electron localization function (ELF) ranges between and 1. In a structural complex, it pro- 
vides a good description of tlic polycentric bonding as a function of the real-space coordinates. 
Generally speaking, a higher ELF implies a lower Pauli kinetic energy, which corresponds to 
having localized covalent bonds or lone electron pairs (ELF = 1 would correspond to a perfect 
localization). For more details, see A. D. Becke and K. E. Edgecombe, J. Chem. Phys. 92, 
5397 (1990); S. Noury, F. Colonna, A. Savin, and B. Silvi, J. Mol. Struct. 450, 59 (1998). 

[24] E. Sanville, S. D. Kenny, R. Smith, and G. Henkelman, J. Comp. Chem. 28, 899 (2007). 



9 



